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ABSTRACT 

We  compare  mid-infrared  emission-line  properties  from  high-resolution  Spitzer  spectra  of  a  hard  X-ray 
(14-195  keV)  selected  sample  of  nearby  (z  <  0.05)  active  galactic  nuclei  (AGNs)  detected  by  the  Burst  Alert 
Telescope  (BAT)  aboard  Swift.  The  luminosity  distribution  for  the  mid-infrared  emission  lines,  [Oiv]  25.89  /xm, 

[Ne  ii]  12.81  /xm,  [Ne  hi]  15.56  /xm,  and  [Ne  v]  14.32/24.32  /xm,  and  hard  X-ray  continuum  show  no  differences 
between  Seyfert  1  and  Seyfert  2  populations;  however,  six  newly  discovered  BAT  AGNs  are  under-luminous  in 
[O  iv],  most  likely  the  result  of  dust  extinction  in  the  host  galaxy.  The  overall  tightness  of  the  mid-infrared  correla¬ 
tions  and  BAT  fluxes  and  luminosities  suggests  that  the  emission  lines  primarily  arise  in  gas  ionized  by  the  AGNs. 

We  also  compare  the  mid-infrared  emission  lines  in  the  BAT  AGNs  with  those  from  published  studies  of  ULIRGs, 
Palomar-Green  quasars,  star-forming  galaxies,  and  LINERs.  We  find  that  the  BAT  AGN  sample  falls  into  a  distinc¬ 
tive  region  when  comparing  the  [Ne  m]/[Ne  ii]  and  the  [O  iv]/[Ne  hi]  ratios.  These  line  ratios  are  lower  in  sources 
that  have  been  previously  classified  in  the  mid-infrared/optical  as  AGNs  than  those  found  for  the  BAT  AGNs,  sug¬ 
gesting  that,  in  our  X-ray  selected  sample,  the  AGNs  represent  the  main  contribution  to  the  observed  line  emission. 

These  ratios  represent  a  new  emission  line  diagnostic  for  distinguishing  between  AGNs  and  star-forming  galaxies. 
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1.  INTRODUCTION 

Active  galactic  nuclei  (AGNs)  span  over  7  orders  of  mag¬ 
nitude  in  bolometric  luminosity  (Lbofi  Koratkar  &  Blaes  1999) 
and  yet  are  all  believed  to  be  powered  by  the  same  physical 
mechanism:  accretion  of  matter  onto  supermassive  black  holes 
(e.g.,  Rees  1984;  Peterson  et  al.  2004).  One  way  to  approach  the 
study  of  AGNs  is  to  concentrate  on  those  in  the  local  universe 
(e.g.,  z  <  0.05),  which  permits  us,  among  other  things,  to  de¬ 
termine  the  properties  of  the  host  galaxy.  Such  studies  tend  to 
focus  on  Seyfert  galaxies,  which  are  modest  luminosity  AGNs 
(Lboi  <  1045  erg  s_1),  but  bright  enough,  due  to  their  proximity, 
to  be  studied  across  the  full  electromagnetic  spectrum. 

Although  Seyfert  galaxies  and  other  AGNs  have  been  tradi¬ 
tionally  defined  in  terms  of  their  optical  properties  (e.g.,  classi¬ 
fication  into  Type  I  and  Type  II;  Khachikian  &  Weedman  1974), 
sample  selection  of  AGNs  via  a  single  waveband  can  lead  to 
observational  bias  (e.g.,  Mulchaey  et  al.  1994).  For  example, 
most  AGNs  are  obscured  from  our  line  of  sight  by  dust  and  gas 
(Matt  2000)  and  any  selection  based  on  optical  (or  UV)  prop¬ 
erties  would  miss  many  objects  or  could  highly  skew  a  sample 
toward  unobscured  objects  (e.g.,  Barger  et  al.  2005).  The  soft 
X-ray  properties  of  Seyfert  galaxies  generally  follow  the  same 
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dichotomy  as  their  optical  properties.  The  X-ray  continuum 
source  in  Seyfert  Is  can  be  observed  directly  (e.g.,  George  et  al. 
1998),  while  the  central  X-ray  source  is  sometimes  undetectable 
in  Seyfert  2s,  due  to  material  with  Ah  >  1022,5  cm-2  along  our 
line  of  sight.  There  is  additionally  a  wide  range  of  effective  IR 
colors  (Kuraszkiewicz  et  al.  2003;  Lutz  et  al.  2005)  which  can 
introduce  selection  bias.  A  comparison  of  infrared  and  X-ray 
data  (Franceschini  et  al.  2005)  shows  a  factor  of  30  range  in 
the  IR  24  /xm  to  ~4  keV  X-ray  flux  ratio  for  X-ray  selected 
AGNs,  suggesting  a  range  of  geometries  and  optical  depths  for 
dust  reprocessing,  and  probably  variance  in  the  intrinsic  power- 
law  AGN  continuum  (see  McKernan  et  al.  2009  for  discussion). 
Obscuration  and  star  formation  in  the  host  galaxy  can  also  dom¬ 
inate  and  introduce  confusion  in  the  IR  (e.g..  Lutz  et  al.  2004; 
Barmby  et  al.  2006).  In  fact,  virtually  all  surveys  for  AGNs 
based  purely  on  IR,  optical,  UV,  or  soft  X-ray  data  have  been 
biased  (Mushotzky  2004).  Even  Sloan  surveys  (Heckman  et  al. 
2004)  or  IR  surveys  (Franceschini  et  al.  2005)  have  required 
indirect  AGN  indicators  which  are  known  to  be  not  necessarily 
robust  (Melendez  et  al.  2008b). 

To  understand  the  intrinsic  properties  of  AGNs  as  a  class, 
it  is  critical  to  start  with  a  survey  where  we  can  be  as  certain 
as  possible  that  we  are  viewing  the  AGN-only  parts  of  these 
galaxies.  At  X-ray  energies  of  E  >  10-20  keV,  the  obscuring 
material  is  relatively  optically  thin  for  column  densities  less 
than  ~3  x  1024  cm-2  (Compton-thin  objects).  Even  if  an  AGN 
is  well  buried  within  its  host  galaxy  there  is  an  unaffected  view 
of  the  central  power  source.  A  hard  X-ray  survey  should  thus 
find  all  Compton  thin  AGNs  in  a  uniform  fashion  and  is  the 
most  representative,  since  at  present,  there  are  very  few,  if  any, 
known  X-ray  “quiet”  AGNs.  Such  a  hard  X-ray  survey  is  now 
available  from  the  Swift  Burst  Alert  Telescope  (BAT).  The  Swift 
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BAT  is  sensitive  over  ~85%  of  the  sky  to  a  flux  threshold  of 
2  x  1CT11  erg  cm-2  s-1  in  the  14-195  keV  band  (Markwardt 
et  al.  2005).  The  BAT  data  are  about  10  times  more  sensitive 
than  the  previous  hard  X-ray  all  sky  survey  (Levine  et  al.  1984). 
The  BAT  detects  all  bright  AGNs,  whether  they  are  obscured  or 
not.  Moreover,  several  of  the  BAT  sources  are  newly  discovered 
AGNs  or  have  been  poorly  studied,  if  at  all,  at  other  wavelengths 
(Winter  et  al.  2008;  Tueller  et  al.  2008,  2010). 

Nevertheless,  although  all  of  the  BAT-detected  objects  are  true 
AGNs,  in  order  to  fully  explore  the  properties  of  these  AGNs 
one  needs  to  take  a  multi-wavelength  approach.  For  example, 
studying  the  IR  properties  of  the  BAT  AGNs  will  provide  insight 
into  the  IR/X-ray  scatter  and  thus  determine  the  true  distribution 
of  IR  properties.  There  have  been  a  large  number  of  studies 
of  the  mid-infrared  emission  line  properties  of  active  galaxies 
using  both  Infrared  Space  Observatory  (Kessler  et  al.  1996) 
and  Spitzer  Space  Telescope  (Werner  et  al.  2004).  The  ratios 
of  high-  and  low-ionization  mid-infrared  emission  lines  have 
been  widely  used  to  separate  the  relative  contribution  of  the 
AGNs  and  star  formation  (e.g.,  Genzel  et  al.  1998;  Sturm  et  al. 
2002;  Dale  et  al.  2006;  Armus  et  al.  2007;  Farrah  et  al.  2007; 
Melendez  et  al.  2008a;  Baum  et  al.  2010).  More  recently,  (Hao 
et  al.  2009,  H09)  used  new  high-resolution  Spitzer  spectroscopy 
to  probe  the  utility  of  mid-infrared  emission  line  diagnostics  as 
a  way  to  separate  active  galaxies  from  star-forming  galaxies. 
In  our  first  study  of  mid-infrared  properties  of  the  BAT  AGNs 
(Melendez  et  al.  2008b),  we  found  the  [Oiv]25.89  // m  to  be  an 
accurate  indicator  of  the  AGN  luminosity,  with  an  uncertainty 
of  ~0.3  dex;  this  result  has  been  confirmed  using  larger  samples 
(Rigby  et  al.  2009;  Diamond-Stanic  et  al.  2009).  Using  a 
complete,  volume-limited,  sample  of  galaxies  (Goulding  & 
Alexander  2009,  GA09)  demonstrated  the  utility  of  high- 
ionization  mid-infrared  emission  lines,  such  as  [Ne  v]  14.32  /xm, 
to  identify  AGNs  including  those  that  were  not  identified  as 
AGNs  in  optical  studies  (see  also,  Landi  et  al.  2007;  Abel  & 
Satyapal  2008;  Dudik  et  al.  2009;  Satyapal  et  al.  2009).  Similar 
results  have  been  found  by  (Bernard-Salas  et  al.  2009,  B09)  in 
their  study  of  starburst  galaxies. 

This  paper  is  the  first  in  a  series  seeking  to  understand  the 
nature  of  the  observed  mid-infrared  luminosities  in  AGNs  and 
their  wide  variety  of  spectral  forms.  Here  we  report  results  from 
the  portion  of  our  sample  that  have  high-resolution  Spitzer  spec¬ 
tra.  This  work  complements  the  extensive  optical  imaging  and 
spectroscopy  of  the  AGN  population  (Koss  et  al.  2010;  Winter 
et  al.  2010)  and  the  detailed  analysis  of  the  X-ray  properties  of 
the  BAT  AGN  sample  (e.g..  Winter  et  al.  2008,  2009a,  2009b). 
In  following  papers,  we  will  report  on  the  results  from  our  anal¬ 
ysis  of  the  low-resolution  Spitzer  spectra  which  will  include 
the  study  of  polycyclic  aromatic  hydrocarbon  features,  silicate 
absorption,  and  mid-infrared  continuum  properties  of  the  BAT 
AGN  sample.  In  order  to  calculate  the  luminosities  presented 
in  this  work  we  assumed  a  flat  universe  with  a  Hubble  constant 
Ha  =  71  km  s-1  Mpc~*,QA  =  0.73,  andQjyj  =  0.27,  withred- 
shift  values  taken  from  NASA’s  ExtraGalactic  Database  (NED), 
except  for  sources  with  redshift  values  of  z  <  0.01  where 
distances  are  taken  from  The  Extragalactic  Distance  Database 
(EDD;  Tully  1988;  Tully  et  al.  2009). 

2.  OBSERVATIONS  AND  ANALYSIS 
2.1.  X-ray  Data  and  AGN  Sample  Selection 

Our  sample  was  selected  from  the  first  unbiased  local  AGN 
sample  obtained  with  the  Swift  BAT  survey.  The  X-ray  positions 


have  ~3"  positional  uncertainties  and  we  are  highly  confident 
of  the  positions  and  the  optical  identifications.  Source  identi¬ 
fications  are  based  primarily  on  the  X-ray  imaging  data  and  a 
correlation  with  optical  images  and  catalogs.  In  some  cases,  the 
identifications  are  based  on  positional  coincidences  with  previ¬ 
ously  known  AGNs. 

The  median  redshift  of  the  BAT  objects  is  z  ~  0.025.  Our 
selection  criteria  are  z  <  0.05,  \b\  >  19°  and  a  hard  X-ray 
BAT  flux  of  >2  x  10~n  erg  cm-2  s  1 .  The  flux  limit  provides 
sufficient  S  /N  to  be  sure  that  the  sources  are  statistically  robust. 
Our  total  BAT  sample  contains  130  objects  above  a  significance 
threshold  of  5.0a.  We  obtained  from  the  Infrared  Spectrograph 
(IRS;  Houck  et  al.  2004)  on  board  Spitzer  high-  and  low- 
resolution  spectra  for  60  of  these  objects,  while  70  previously 
had  Spitzer  IRS  spectra  from  other  observing  programs.  There 
are  published  BAT  AGN  catalogs  from  two  surveys,  the  9  month 
survey  (Tueller  et  al.  2008)  and  an  updated  22  month  catalog 
(Tueller  et  al.  2010).  Here  we  use  the  22  month  fluxes  from 
Tueller  et  al.  (2010)  as  the  22  month  survey  is  the  most  sensitive. 
However,  in  our  sample,  four  of  the  AGNs  in  the  9  month  survey 
were  not  detected  in  the  22  month  survey;  therefore,  we  will  use 
the  fluxes  from  the  9  month  survey  as  upper  limits  (see  footnote 
in  Table  1). 

2.2.  Spitzer  Observations  and  Data  Analysis 

In  this  work,  we  discuss  only  the  high-resolution  Spitzer 
spectroscopy  of  our  total  sample,  which  result  on  a  subsample 
of  79  BAT  AGNs.  The  galaxies  in  our  sample  were  observed 
with  IRS  in  the  short-high  (SH,  X  =  9.9-19.6  pm,  4"7  x  1 1"3, 
R  ~  600)  and  long-high  (LH,  X  =  18.7-37.2  pm,  1 1"1  x  22"3, 
R  ~  600)  IRS  order  in  the  staring  mode.  Science  targets  in  the 
staring  mode  are  placed  at  two  node  position  along  the  IRS  slit. 
This  sample  includes  fluxes  found  in  the  literature  (Weedman 
et  al.  2005;  Armus  et  al.  2006;  Tommasin  et  al.  2008)  and  from 
our  analysis  of  unpublished  archival  spectra  observed  with  IRS, 
including  observations  collected  with  the  Spitzer  Cycle  3  and 
Cycle  5  GTO  program  “ Spitzer  Observations  to  Complete  the 
First  Unbiased  AGN  Sample  of  the  Local  Universe”  (P30745 
and  P50588,  PI:  K.  Weaver).  This  high-resolution  sample 
includes  38  Seyfert  Is,  33  Seyfert  2s,  six  previously  unknown 
or  poorly  studied  AGNs  and  two  LINERs.  High-resolution 
spectroscopy  allows  us  to  easily  separate  blended  features 
such  as  [Oiv]  — [Fen]  25.99 pm  and  [Nev]  14.32 pm—  [Cl n] 
14.37  pm.  Moreover,  the  subsample  presented  in  this  work 
spans  the  same  range  on  14-195  keV  luminosities  as  the  22 
month  BAT  AGN  sample  and  is  large  enough  to  be  statistically 
representative  of  the  whole  sample  of  130  AGNs.  Furthermore, 
the  Kolmogorov-Smirnov  (K-S)  test  probabilities  for  Seyfert 
1  and  Seyfert  2  galaxies  between  the  subsample  presented  in 
this  work  and  the  whole  sample  of  130  sources  are  46.2%  and 
11. 5%, 13  respectively,  meaning  that  there  are  no  differences  in 
their  X-ray  luminosities.  Finally,  our  sample  also  has  the  same 
range  in  mid-infrared  emission  line  luminosities  as  complete 
samples  such  as  the  12  pm  (Tommasin  et  al.  2010)  and  the 
revised  Shapley-Ames  sample  (Rigby  et  al.  2009). 

For  the  analysis  of  Spitzer  data  we  used  the  basic  calibrated 
data  files  preprocessed  using  the  S17.2  IRS  pipeline.  This 
includes  ramp  fitting,  dark  sky  subtraction,  drop  correction, 

13  A  probability  value  of  less  than  5%  represents  a  high  level  of  significance 
that  two  samples  drawn  from  the  same  parent  population  would  differ  this 
much  5%  of  the  time,  i.e.,  that  they  are  different.  A  strong  level  of  significance 
is  obtained  for  values  smaller  than  1%  (e.g.,  Press  et  al.  1992;  Bevington  & 
Robinson  2003). 
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Table  1 

Spitzer  IRS  High-resolution  Spectroscopy  of  the  BAT  Sample  of  AGNs 


Name 

Distance  Type3 
(Mpc) 

BATb 

[Neil] 

(12.81  /xm) 

[Nev]  [Nem]  [Nev] 

(14.32 /xm)  (15.56 /xm)  (24.32 /xm) 

[Oiv] 

(25.89  (tm) 

Integrated  Line  Fluxes  (10-21  W  cm-2) 

2MASX  J05580206-3820043 

146.8 

1 

3.99  ±  0.37 

3.20  ±  0.89 

1.89  ±  0.13 

4.72  ±  0.14 

3.80  ±  0.51 

3.69  ±  0.45 

3C120 

143.0 

1 

11.89  ±  0.63 

9.19  ±  0.66 

17.21  ±  0.28 

27.43  ±  0.92 

27.63  ±  8.92 

116.68  ±  1.22 

Ark  120 

141.7 

1 

7.08  ±  0.57 

3.47  ±  0.91 

<1.18 

4.33  ±  0.79 

<1.29 

4.03  ±  0.38 

Cen  Ac-d 

4.9 

2 

92.62  ±  0.71 

193.00 

23.16 

147.65 

29.92 

131.24 

ESO  005-G004d 

22.4 

4.48  ±  0.60 

16.64  ±  1.19 

0.50  ±  0.12 

5.29  ±  0.34 

<0.49 

4.49  ±  0.19 

ESO  033-  G002 

77.5 

2 

2.61  ±  0.45 

2.67  ±  0.18 

5.57  ±  0.56 

8.34  ±  0.56 

5.37  ±  0.51 

13.85  ±  0.73 

ESO  103-G035 

56.7 

2 

1.1.14  :  0.59 

30.92  ±  2.07 

15.85  ±  3.49 

41.62  ±  1.78 

10.31  ±  0.87 

34.40  ±  0.75 

ESO  140-G043 

60.5 

1 

4.57  ±  0.66 

11.41  ±  0.26 

7.83  ±  0.70 

13.93  ±  0.32 

7.95  ±  0.44 

27.46  ±  0.35 

ESO  323-G077 

64.2 

1.2 

4.70  ±  0.66 

40.16  ±  1.41 

5.44  ±  1.61 

18.07  ±  1.07 

6.19  ±  0.14 

25.11  ±  0.62 

ESO  362-GO  18 

53.1 

1.5 

6.22  ±  0.52 

9.94  ±  0.23 

3.15  ±  0.37 

7.16  ±  0.18 

2.75  ±  0.20 

10.11  ±  0.37 

ESO  417-  G006 

69.7 

2 

3.06  ±  0.46 

1.77  ±  0.14 

0.42  ±  0.10 

2.54  ±  0.07 

0.47  ±  0.16 

4.04  ±0.18 

F9 

205.9 

1 

5.07  ±  0.45 

2.91  ±  0.68 

<2.66 

5.31  ±  0.37 

2.74  ±  0.37 

6.08  ±  0.26 

F49 

86.2 

2 

2.93  ±  0.54 

39.55  ±  2.93 

25.50  ±  2.96 

47.12  ±  2.43 

12.74  ±  0.36 

39.22  ±  2.37 

IC  486 

114.8 

1 

3.22  ±  0.70 

6.86  ±  0.05 

3.30  ±  0.20 

6.79  ±  0.21 

2.84  ±  0.01 

11.19  ±  0.29 

IC  1816 

72.5 

1 

2.58  ±  0.48 

18.36  ±  1.51 

6.21  ±  0.61 

19.64  ±  1.17 

4.34  ±  0.30 

16.52  ±  1.81 

IC  4329A 

68.7 

1.2 

33.08  ±  0.62 

30.19  ±  2.76 

32.79  ±  6.01 

65.30  ±  0.95 

29.90  ±  1.74 

103.58  ±  3.24 

IC  5063 

48.3 

2 

8.59  ±  0.72 

28.22  ±  3.34 

35.97  ±  1.58 

73.67  ±  4.61 

25.86  ±  1.20 

117.21  ±  11.13 

MCG-01-13-025e 

68.0 

1.2 

<4.5 

2.87  ±  0.18 

<0.18 

1.90  ±  0.05 

<1.44 

0.76  ±0.18 

MCG-05-23-016 

36.1 

2 

20.77  ±  0.56 

18.13  ±  0.42 

11.06  ±  0.50 

16.96  ±  0.93 

11.19  ±  0.26 

27.95  ±  7.61 

MCG-01-24-012 

84.2 

2 

4.58  ±  0.51 

6.74  ±  0.41 

2.90  ±  0.27 

6.07  ±  0.36 

2.30  ±  0.58 

10.07  ±  0.77 

MCG-02-58-22 

205.1 

1.5 

10.17  ±  0.57 

5.29  ±  0.39 

2.80  ±  0.34 

9.71  ±  0.16 

3.05  ±  0.39 

13.74  ±  1.74 

MCG-03-34-064 

70.8 

1.8 

3.15  ±  0.45 

57.83  ±  7.29 

54.98  ±  5.26 

120.13  ±  6.88 

35.46  ±  1.92 

110.83  ±  9.35 

MCG-06-30-015 

32.9 

1.2 

7.82  ±  0.57 

4.20  ±  0.12 

6.05  ±  0.76 

6.62  ±  0.33 

6.62  ±  0.05 

23.49  ±  0.81 

MRK  3C 

57.7 

2 

15.65  ±  0.61 

100.00 

64.50 

179.00 

67.50 

214.00 

MRK  6f 

80.6 

1.5 

7.61  ±  0.55 

28.00  ±  0.23 

9.39  ±  0.19 

49.34  ±  0.32 

10.43  ±  0.21 

48.24  ±  0.27 

MRK  18e 

47.2 

<3.1 

17.04  ±  0.71 

0.78  ±  0.01 

8.24  ±  0.24 

<2.25 

2.69  ±  0.14 

MRK  79 

95.3 

1.2 

4.89  ±  0.52 

11.29  ±  3.71 

9.62  ±  1.07 

20.43  ±  0.50 

13.16  ±  1.10 

39.99  ±  2.68 

MRK  335 

111.1 

1.2 

2.47  ±  0.41 

2.10  ±  0.31 

1.31  ±  0.11 

2.70  ±  0.15 

<1.55 

6.31  ±  0.17 

MRK  348 

64.2 

2 

13.66  ±  0.56 

15.34  ±  0.74 

6.90  ±  0.36 

20.60  ±  0.79 

4.75  ±  0.27 

17.87  ±  0.23 

MRK  352 

63.5 

1 

4.16  ±  0.50 

0.54  ±  0.08 

<0.16 

0.34  ±  0.02 

<0.23 

0.26  ±  0.04 

MRK  509 

149.2 

1.2 

9.44  ±  0.68 

11.98  ±  1.06 

7.61  ±  1.74 

17.32  ±  3.29 

7.60  ±  0.10 

27.54  ±  0.42 

MRK  590c 

112.7 

1.2 

<3.7 

5.26  ±  0.46 

<3.22 

5.42  ±  0.12 

<1.23 

1.79  ±  0.26 

MRK  766 

55.2 

1.5 

2.42  ±  0.29 

24.27  ±  0.94 

22.14  ±  0.23 

23.52  ±  1.41 

18.32  ±  1.60 

46.47  ±  0.84 

MRK  817 

136.1 

1.5 

2.21  ±  0.36 

4.57  ±  0.88 

2.88  ±  0.48 

5.51  ±  0.77 

4.22  ±  0.73 

6.06  ±  0.23 

MRK  841 

158.2 

1.5 

2.93  ±  0.37 

3.18  ±  0.18 

8.12  ±  0.31 

12.46  ±  0.73 

6.96  ±  0.33 

22.91  ±  1.08 

NGC  454e 

51.8 

2 

<2.3 

4.70  ±  0.40 

3.73  ±  0.14 

6.40  ±  0.17 

5.54  ±  0.31 

15.80  ±  0.85 

NGC513 

83.8 

2 

2.06  ±  0.43 

10.38  ±  1.78 

1.52  ±  0.23 

4.71  ±  0.32 

1.38  ±  0.21 

6.69  ±  0.41 

NGC  788 

58.1 

2 

9.33  ±  0.57 

6.59  ±  0.11 

5.35  ±  0.34 

13.60  ±  0.50 

8.68  ±  1.28 

24.14  ±  0.36 

NGC  931 

71.2 

1.5 

6.56  ±  0.50 

5.50  ±  0.59 

9.91  ±  2.31 

15.94  ±  0.52 

13.52  ±  0.58 

42.90  ±  0.63 

NGC  973 

69.3 

3.09  ±  0.58 

8.03  ±  0.09 

2.07  ±  0.18 

8.92  ±  0.38 

3.21  ±  0.27 

14.11  ±  2.64 

NGC  1052d 

17.8 

LINER 

3.75  ±  0.67 

22.22  ±  0.94 

0.64  ±  0.03 

12.31  ±  0.48 

2.69  ±  1.33 

2.37  ±  0.10 

NGC  1194 

58.0 

1.9 

3.64  ±  0.60 

3.27  ±  0.34 

4.48  ±  0.12 

8.27  ±  0.41 

3.98  ±  0.32 

14.39  ±  0.10 

NGC  1365d 

18.9 

1.8 

7.19  ±  0.44 

161.67  ±  17.48 

22.35  ±  1.97 

61.06  ±  0.90 

38.53  ±  1.50 

145.38  ±  8.89 

NGC  2110 

33.1 

2 

35.01  ±  0.70 

60.19  ±  5.34 

5.22  ±  0.82 

47.40  ±  0.71 

7.65  ±  0.63 

45.71  ±  3.41 

NGC  2992d 

30.5 

2 

4.82  ±  0.63 

53.65  ±  3.66 

32.62  ±  5.38 

61.06  ±  1.98 

27.81  ±  0.20 

114.22  ±  6.44 

NGC  3079c,d 

20.4 

2 

3.44  ±  0.44 

104.00 

1.04 

22.88 

<1.56 

9.26 

NGC  308  ld 

32.5 

2 

10.24  ±  0.67 

12.62  ±  1.16 

12.62  ±  0.81 

36.46  ±  1.25 

35.79  ±  0.05 

119.73  ±  8.42 

NGC  3227d 

20.6 

1.5 

14.13  ±  0.50 

65.05  ±  6.91 

26.23  ±  1.23 

74.62  ±  2.00 

17.82  ±  0.99 

64.91  ±  2.31 

NGC  3281 

45.5 

2 

9.01  ±  0.66 

19.94  ±  2.16 

47.51  ±  2.75 

58.35  ±  2.48 

42.15  ±  3.24 

174.65  ±  13.25 

NGC  3516d,f 

38.9 

1.5 

12.54  ±  0.45 

8.07  ±  0.25 

7.88  ±  0.50 

17.72  ±  0.33 

10.39  ±  0.33 

46.92  ±  0.35 

NGC  3783d 

38.5 

1 

19.45  ±  0.66 

19.82  ±  0.79 

19.54  ±  3.44 

26.18  ±  0.57 

13.82  ±  0.29 

39.25  ±  0.07 

NGC  405  ld 

17.1 

1.5 

4.34  ±  0.35 

19.69  ±  0.92 

10.17  ±  0.65 

16.35  ±  0.44 

16.19  ±  2.53 

36.95  ±  2.05 

NGC  4102d 

17.0 

LINER 

2.58  ±  0.38 

349.81  ±  39.88 

<8.10 

30.44  ±  1.12 

<6.29 

20.66  ±  7.69 

NGC  4138d 

17.0 

1.9 

3.69  ±  0.45 

3.23  ±  0.52 

0.54  ±  0.10 

2.82  ±  0.37 

<0.83 

2.03  ±0.18 

NGC  415  lc>d 

20.3 

1.5 

62.23  ±  0.46 

134.00 

77.72 

204.35 

67.67 

236.51 

NGC  4235d 

35.1 

1 

2.40  ±  0.55 

3.69  ±  0.65 

<0.58 

3.70  ±  0.54 

0.92  ±  0.19 

3.77  ±  0.77 

NGC  4388d 

16.8 

2 

34.64  ±  0.52 

79.74  ±  4.76 

45.35  ±  0.84 

108.18  ±  1.56 

64.10  ±  0.12 

311.42  ±  25.79 

NGC  4395d 

3.9 

1.8 

3.12  ±  0.41 

4.74  ±  0.25 

1.47  ±  0.13 

7.02  ±  0.58 

1.49  ±  0.62 

8.16  ±  0.26 

NGC  4507 

50.3 

2 

22.51  ±  0.68 

33.73  ±  2.63 

12.50  ±  1.39 

28.63  ±  2.36 

10.72  ±  1.96 

36.33  ±  4.14 

NGC  4593d 

39.5 

1 

9.79  ±  0.62 

8.31  ±  0.30 

5.56  ±  1.27 

7.89  ±  0.59 

8.02  ±  1.96 

13.39  ±  2.00 

NGC  4686 

71.6 

3.08  ±  0.45 

2.13  ±  0.16 

<0.22 

1.71  ±  0.13 

0.53  ±  0.09 

0.91  ±  0.30 

NGC  526A 

81.9 

1.5 

5.96  ±  0.51 

6.22  ±  1.61 

3.94  ±  0.47 

9.59  ±  1.10 

3.87  ±  0.86 

17.34  ±  0.84 

NGC  5506d 

28.7 

1.9 

25.64  ±  0.50 

91.75  ±  3.31 

58.28  ±  3.34 

152.13  ±  9.13 

63.25  ±  2.40 

252.82  ±  2.29 

NGC  5548 

73.5 

1.5 

8.08  ±  0.50 

8.93  ±  0.70 

1.95  ±  0.10 

8.99  ±  0.97 

1.73  ±  0.75 

12.75  ±  0.80 

1154 
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Table  1 

(Continued) 

Name  Distance  Typea  BATb  [Neil]  [Nev]  [Nem]  [Nev]  [Oiv] 

(Mpc)  (12.81  /xm)  (14.32  (tm)  (15.56  /tm)  (24.32 /xm)  (25.89  /xm) 


Integrated  Line  Fluxes  (10  21  W  cm  2) 


NGC  5728d 

42.2 

2 

10.54 

± 

0.71 

30.44 

± 

1.81 

23.56 

± 

0.77 

54.76 

± 

0.51 

28.15  ±  0.56 

118.40 

± 

7.61 

NGC  5995 

108.5 

2 

4.51 

± 

0.61 

13.32 

± 

1.91 

7.60 

± 

1.42 

9.19 

± 

0.59 

4.61  ±  0.84 

12.20 

d= 

0.80 

NGC  6240E 

105.4 

2 

7.30 

± 

0.62 

193.10 

± 

3.70 

5.10 

± 

0.90 

70.40 

± 

2.40 

<3.90 

27.20 

± 

0.70 

NGC  6860 

63.6 

1 

6.50 

± 

0.73 

5.88 

± 

0.11 

2.95 

± 

0.94 

7.50 

± 

0.53 

2.84  ±  0.21 

11.70 

± 

0.39 

NGC  7172d 

33.9 

2 

18.11 

± 

0.70 

32.03 

± 

2.44 

8.99 

± 

0.80 

15.71 

± 

0.69 

11.26  ±  0.95 

42.60 

± 

3.44 

NGC  7213d 

22.0 

1.5 

5.75 

± 

0.67 

27.47 

± 

1.33 

0.56 

± 

0.01 

12.77 

± 

0.66 

<1.94 

2.75 

± 

0.59 

NGC  7314d 

18.3 

1.9 

4.63 

± 

0.59 

8.97 

± 

0.74 

16.84 

± 

0.60 

23.28 

± 

0.41 

21.92  ±  0.40 

69.62 

± 

7.82 

NGC  7469c 

69.8 

1.2 

6.66 

± 

0.44 

200.00 

11 

1.60 

34 

.00 

14.70 

43.00 

NGC  7582d 

22.0 

2 

7.92 

± 

0.55 

250.94 

± 

3.53 

38.02 

± 

3.86 

104.99 

± 

4.37 

59.54  ±  6.39 

227.65 

± 

10.88 

NGC  7603 

127.6 

1.5 

4.70 

± 

0.51 

9.32 

± 

0.72 

0.47 

± 

0.01 

5.07 

± 

0.49 

<0.38 

3.34 

± 

0.18 

NGC  7682 

73.4 

2 

2 

.27 

5.46 

± 

0.25 

1.98 

± 

0.19 

8.07 

± 

0.15 

3.69  ±  1.46 

16.21 

± 

0.66 

UGC  03601 

73.3 

1.5 

4.38 

± 

0.67 

5.26 

± 

0.53 

2.01 

± 

0.12 

7.60 

± 

0.13 

1.79  ±  0.44 

13.45 

± 

0.65 

UGC  06728 

27.7 

1.2 

2.95 

± 

0.37 

1.40 

± 

0.36 

1.16 

± 

0.10 

1.38 

± 

0.05 

3.24  ±  0.03 

4.61 

± 

0.81 

UGC  12282 

72.7 

2.49 

± 

0.50 

2.08 

± 

0.28 

0.50 

± 

0.01 

1.76 

it 

0.05 

1.44  d=  0.20 

4.07 

± 

0.37 

UGC  12741 

74.7 

4.00 

± 

0.59 

1.72 

± 

0.07 

0.39 

± 

0.06 

1.55 

± 

0.05 

0.65  ±  0.06 

3.58 

± 

0.38 

Notes. 

a  AGN  types  are  taken  from  NED. 

b  The  BAT  flux  is  presented  in  units  of  10  11  erg  s-1  cm~2. 
c  Weedman  et  al.  (2005). 
d  Distances  are  taken  from  EDD. 
c  From  the  9  month  BAT  survey  (Tueller  et  al.  2008). 
f  Tommasin  et  al.  (2008). 
g  Armus  et  al.  (2006). 


linearity  correction,  and  wavelength  and  flux  calibrations.14 
Many  of  the  sources  in  our  sample  have  dedicated  off-source 
observations  to  do  sky  subtraction,  which  can  alleviate  the  effect 
of  rogue  pixels  and  variable  background.  Off-source  images 
were  averaged  for  each  node  position  and  order  to  obtain  a 
final  background  image.  Then,  we  subtract  the  sky  background 
from  our  spectrum.  The  full  slit  spectra  were  extracted  from 
the  IRS  data  using  the  Spectroscopy  Modeling  Analysis  and 
Reduction  Tool  (SMART)  v6.4.0  (Higdon  et  al.  2004).  For 
the  extraction  we  used  the  “Full  Aperture”  extraction  method 
for  high-resolution  observations  of  point  sources.  We  created 
median  basic  calibrated  data  files  from  each  node  and  then 
the  spectra  from  each  node  position  for  SH/LH  were  averaged 
using  2.5a  clipping  to  reject  outliers.  Finally,  we  trimmed  the 
edges  of  the  orders  to  obtain  a  clean  spectrum.  For  sources 
without  dedicated  off-source  observations  we  did  not  perform 
any  background  subtraction  because  we  only  required  emission¬ 
line  fluxes;  furthermore,  our  hard  X-ray  selected  sample  is 
characterized  by  bright  nuclear  sources  which  fill  the  high- 
resolution  slit  resulting  in  a  minimal  background  correction. 
We  performed  the  line  fit  with  SMART  using  a  polynomial 
to  fit  the  continuum  and  a  Gaussian  for  the  line  profile.  In 
Table  1,  we  report  the  line  fluxes,  together  with  their  la 
statistical  error  for  the  whole  high-resolution  sample.  The 
typical  la  statistical  errors  for  the  mission  lines  presented  in 
this  work  are,  on  average,  ~9%,  ~10%,  ~5%,  ~12%,  and  ~8% 
for  the  [Ne  ii]  12.81 /xm,  [Ne  v]  14.32 /xm,  [Ne  hi]  15.56  /im, 
[Nev]  24.32  /xm,  and  [Oiv]  emission  lines,  respectively.  From 
this  and  for  the  sake  of  simplicity  we  did  not  plot  error  bars 
on  individual  objects  in  the  different  comparisons  presented  in 
this  work,  as  the  errors  are  comparable  to  the  symbol  size  uses 
in  the  figures.  For  non-detections  we  quote  the  3a  upper  limits 


14  See  the  IRS  Pipeline  Handbook,  http://ssc.spitzer.caltech.edu/irs/. 


as  defined  for  emission  lines  with  an  S/N  less  than  3.  Finally, 
the  emission  line  fluxes  are  presented  here  without  reddening 
corrections. 

Each  of  the  galaxies  in  our  sample  has  clearly  detected 
[O  iv],  [Ne  ii],  and  [Ne  hi]  emission.  However,  [Ne  v]  14.32  /xm 
and  24.32 /xm  were  not  detected  in  ~10%  and  ~15%  of 
the  sources,  respectively.  Similar  results  are  obtained  for  the 
[Ne  v]  14.32  /xm  and  24.32  /xm  in  the  Spitzer  high-resolution 
spectroscopy  of  the  12  /xm  Seyfert  galaxies  (Tommasin  et  al. 
2008,  2010).  The  importance  of  [Oiv]  and  the  neon  emission 
lines  is  that  they  are  sufficient  to  distinguish  between  stellar 
and  AGN  activity.  For  example,  the  [Ne m]/[Ne ii]  and  [Oiv]/ 
[Ne  ii]  ratios  are  good  discriminators  of  star  formation  and  AGN 
emission  (Genzel  et  al.  1998;  Sturm  et  al.  2002;  Satyapal  et  al. 
2004;  Melendez  et  al.  2008a),  while  the  strengths  of  [Oiv] 
(Melendez  et  al.  2008b;  Rigby  et  al.  2009)  and  [Nev]  (e.g., 
Satyapal  et  al.  2007;  Dudik  et  al.  2007;  Abel  &  Satyapal  2008) 
scale  with  the  luminosity  of  the  AGNs.  Moreover,  the  neon  line 
ratios  are  also  tracers  that  are  insensitive  to  abundance.  It  should 
be  noted,  however,  that  [Oiv]  and  [Nev]  can  be  weak  in  some 
classical  AGNs  (Weedman  et  al.  2005),  notably  ultraluminous 
infrared  galaxies  (UFIRGs),  possibly  due  to  the  shielding  of 
the  narrow  line  region  by  optically  thick,  dusty  gas  close  to  the 
AGNs  (e.g.,  Armus  et  al.  2007). 

3.  MID-INFRARED  PROPERTIES  OF  THE  BAT  SAMPFE 
3.1.  Comparison  of  the  IR  Emission  Line  and  BAT  Luminosities 

As  noted  above,  via  the  hard  X-ray  band  we  can  detect 
AGNs  that  may  have  been  missed  in  optical  surveys  due  to  the 
extinction  of  their  optical  emission  lines.  Among  the  79  targets 
presented  here  (BAT  AGNs),  we  have  found  six  poorly  studied  or 
previously  unknown  AGNs  (“new  BAT-detected  AGN”):  ESO 
005-G004,  MRK  18,  NGC  973,  NGC  4686,  UGC  12282,  and 
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Figure  1.  Combined  Spitzer  IRS  SH  and  LH  spectra  for  the  newly  detected  AGNs  in  the  BAT  sample.  The  [Ne  n]  12.81  //m.  [Nem]  15.56 /am.  [Ne  v]  14.32 /am, 
[Ne  v]  24.32  /am,  and  [O  iv]  25.89  /am  lines  are  indicated  by  the  vertical  lines.  The  strong  emission  line  features  at  ~33  /am  and  ~35  /am  correspond  to  [S  m]  and 
[Si  n],  respectively. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  2.  Histograms  for  the  mid-infrared  and  BAT  (14-195  keV)  luminosities  for  our  sample.  Seyfert  1  and  Seyfert  2  galaxies  are  indicated  with  gray  bins  and 
solid  lines,  respectively.  The  newly  detected  BAT  AGNs  are  indicated  in  bins  with  dashed  lines;  note  that  their  [Nem]  and  [Oiv]  luminosities  appear  to  be  low, 
i.e.,  compared  to  AGNs  with  similar  BAT  luminosities.  The  K-S  test  for  these  emission  line  luminosities  shows  that  two  samples  drawn  from  the  same  population 
would  differ  this  much  ~6.3%,  ~82.6%,  ~69.8%,  ~50.9%,  ~28.6%,  and  ~98.1%  of  the  time  for  the  BAT,  [Neil],  [Ne  v]  14  /xm,  [Ne  hi],  [Ne  v]  24  /xm,  and  [Oiv] 
luminosity  distributions,  respectively. 
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Figure  3.  Comparison  between  the  ratios  of  the  BAT  band  and  the  mid-infrared 
emission  lines.  Seyfert  1  galaxies  are  presented  as  black  circles,  Seyfert  2 
galaxies  are  red  triangles,  blue  squares  represent  the  newly  detected  BAT  AGNs 
and  green  stars  are  LINERs.  Of  particular  interest  is  the  statistically  significant 
(see  Table  3)  separation  between  Seyfert  types  in  the  E14_195  key/L[Oiv]  ratio, 
in  agreement  with  this  ratio  being  an  indicator  of  Compton  scattering  in  the 
BAT  band.  One  must  note  that  below  Lbat/[Oiv]  <~  2.0  there  are  only  two 
Seyfert  1  galaxies  but  13  Seyfert  2  galaxies  (~40%  of  the  Seyfert  2  population). 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

UGC  12741.  There  have  been  attempts  to  classify  some  of 
these  galaxies.  Of  these,  ESO  005-G004  has  been  classified 
as  a  Seyfert  2,  based  on  its  X-ray  properties  although  no  optical 
AGN  signature  has  been  detected  (Landi  et  al.  2007;  Ueda  et  al. 

2007) .  Based  on  follow-up  observations  in  the  optical,  MRK 
18  has  an  ambiguous  classification  with  Hii/LINER  properties 
(Winter  et  al.  2010)  and  NGC  973  has  been  classified  as  a  narrow 
emission  line  AGN  with  Sy2/LINER  properties  (Masetti  et  al. 

2008) .  UGC  12282  has  been  classified  as  an  Sy  1.9  based  on 
its  optical  spectra  (Veron-Cetty  &  Veron  2006).  Finally,  to  the 
best  of  our  knowledge,  there  is  no  assigned  classification  for 
NGC  4686  and  UGC  12741.  The  Spitzer  IRS  spectra  of  these 
objects  are  shown  in  Figure  1 . 

In  Figures  2-7,  we  compare  the  mid-infrared  emission  lines 
and  the  BAT  luminosities  (and  fluxes),  correlations  between 
the  emission  lines,  and  emission  line  ratios.  The  statistical 
analysis  for  these  plots  is  listed  in  Table  2,  which  includes  the 
Spearman  rank  order  correlation  coefficient  with  its  associated 
null  probability,  the  generalized  Kendall’s  correlation  coefficient 
for  censored  data  and  the  Kendall’s  coefficient  for  partial 
correlation  with  censored  data.  One  should  note  that,  due  to 
redshift  effects,  luminosity-luminosity  plots  will  almost  always 


show  some  correlation.  Thus,  we  are  primarily  interested  in  the 
dispersion  of  the  correlations  or  the  slopes.  Furthermore,  caution 
must  be  taken  when  applying  statistical  analysis  to  data  sets 
that  contain  non-detections  (upper  limits),  or  “censored”  data 
points.  To  deal  with  these  problems  we  have  used  Astronomy 
SURVival  analysis  (ASURV)  Rev  1 .2  (Isobe  &  Feigelson  1990), 
which  implement  the  methods  presented  in  Isobe  et  al.  (1986). 
We  also  used  a  test  for  partial  correlation  with  censored  data 
(Akritas  &  Siebert  1996)  in  order  to  exclude  the  redshift  effect 
in  the  correlations. 

In  Figure  2,  we  show  the  distribution  of  the  observed  [Ne  ii], 
[Nem],  [Nev],  [Oiv],  and  BAT  luminosities  (That)  for  our 
sample.  The  results  of  the  K-S  tests,  comparing  Seyfert  Is 
and  2s,  for  these  quantities  are  listed  in  Table  3.  This  table 
also  includes  information  about  the  numbers  of  Seyfert  1  and 
Seyfert  2  galaxies,  median  values,  and  standard  deviations  of 
the  mean  for  the  measured  quantities.  As  is  apparent  from 
both  the  histograms  and  the  K-S  results,  the  two  Seyfert  types 
are  statistically  indistinguishable,  i.e.,  there  are  essentially  no 
differences  between  Seyfert  Is  and  2s  when  we  directly  compare 
the  BAT  luminosities  or  any  of  the  IR  line  strengths.  These 
results  are  in  agreement  with  the  similar  optical  luminosities 
found  for  broad  and  narrow  line  sources,  with  optical  Seyfert 
classifications,  in  the  most  recent  study  on  the  optical  spectral 
properties  of  the  BAT  AGNs  (Winter  et  al.  2010).  The  smallest 
null-test  probability  is  for  /.  bat.  which,  as  we  will  discuss  below, 
is  evidence  that  Seyfert  2s  suffer  more  from  Compton  scattering 
in  the  BAT  band.  Although  there  were  too  few  of  the  “new 
BAT  AGN”  to  separate  them  for  a  K-S  on  these  quantities, 
the  histograms  indicate  that  they  are  relatively  weaker  in  their 
observed  mid-infrared  emission  line  luminosities,  compared  to 
their  BAT  luminosities. 

As  mentioned  before,  a  detailed  statistical  analysis  for  the 
different  correlations  between  the  BAT  and  mid-infrared  lumi¬ 
nosities  is  presented  in  Table  2.  From  this  analysis  the  weaker 
correlation  found  in  the  [Ne  ii]-BAT  relationship  could  be  the 
result  of  active  star  formation  contributing  to  the  [Ne  ii]  emission 
line  in  some  of  the  AGNs  (e.g.,  Sturm  et  al.  2002;  Schweitzer 
et  al.  2006;  Ho  &  Keto  2007;  Melendez  et  al.  2008a).  On  the 
other  hand,  stronger  correlations  are  found  between  [Nem], 
[Nev],  and  [Oiv]  luminosities  when  compared  to  /.bat,  sug¬ 
gesting  that,  on  average,  there  is  no  strong  enhancement  due 
to  star  formation  in  the  [Nem]  and  [Oiv]  emission  in  the  BAT 
sample. 

In  Table  4,  we  present  the  linear  regression  fits  for  all 
correlations.  The  somewhat  steeper  slopes  for  the  Seyfert 
2s  are  consistent  with  the  decrement  in  the  observed  X-ray 
emission  in  the  BAT  band  due  to  the  effect  of  Compton  scatter. 
The  strong  similarities  between  Seyfert  types  lessen  when 
we  compare  the  ratios  of  mid-infrared  line  strengths  and  the 
X-ray  continuum  strength.  In  Figure  3,  we  show  the  ratios 
of  Lbat  and  [Nen],  [Nem],  and  [Nev]  14.32 pm  compared 
to  /-bat/[0  iv].  The  histograms  for  the  ratios  of  Lbat  and 
[Ne  ii],  [Ne  m],  [Ne  v]  14.32/24.32  p. m  and  [O  iv]  are  presented 
in  Figure  4.  Of  particular  interest  in  these  comparisons  is 
that  below  log  Lbat/[0  iv]  <~  2.0  there  are  only  two  Seyfert 

1  galaxies  but  13  Seyfert  2  galaxies  (~40%  of  the  Seyfert 

2  population  in  our  sample).  Therefore,  in  these  plots,  one 
can  see  that  the  Seyfert  2s  have  lower  ratios  of  Lbat  to 
the  emission  lines  compared  to  the  Seyfert  Is,  which  is,  we 
believe  is  due  to  the  effect  of  Compton  scattering  in  the  BAT 
band  in  some  Seyfert  2  galaxies.  This  result  is  in  agreement 
with  the  K-S  test  null  probability  for  this  ratio  (see  Table  3). 
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Figure  4.  Histograms  for  the  ratios  of  the  BAT  band  and  the  mid-infrared  emission  lines.  Seyfert  1  and  Seyfert  2  galaxies  are  indicated  with  gray  bins  and  solid 
lines,  respectively.  The  newly  detected  BAT  AGNs  are  indicated  in  bins  with  dashed  lines.  Interestingly,  all  these  ratios  are  log  normal  distributed  with  a  standard 
deviation  of  only  0.5  dex,  except  for  the  L 14-195  kev/^[Nem]  with  a  standard  deviation  of  0.4  dex  (see  Table  3).  The  K-S  test  for  these  ratios  shows  that  two  samples 
drawn  from  the  same  population  would  differ  this  much  ~2.8%,  ~9.0%,  ~2.4%,  ~26.1%  and  ~9.5%  of  the  time  for  the  L 14-195 kev/^[Nen]»  ^  14-195  kev/£[Nem]> 
L 14-195 kev/£[0 iv]>  ^14-195 kev/^[Nev](  14.32 /im),  and  L 14-195 kev/^[Ne v](24.32/im)  distributions,  respectively. 


Similarly,  there  is  only  one  Seyfert  2  galaxy,  NGC  4138,  with 
logLBAT/[Oiv]  >  3.0. 

From  the  Tbat/[Oiv]  ratios,  the  lowest  values,  likely  due 
to  the  effects  of  Compton  scattering  in  the  BAT  band  are 
NGC  1365,  NGC  2992,  NGC  3281,  NGC  7582,  MRK  3, 
and  MCG-03-34-064.  Note  that  NGC  2992  and  NGC  1365 
are  very  variable  in  the  X-ray,  which  can  contribute  to  their 
low  Lbat/[Oiv]  and  TBAr/[Nev]  ratios  (see  Gilli  et  al.  2000; 
Risaliti  et  al.  2000b).  In  the  BAT  spectra  (R.  F.  Mushotzky  et  al. 
2010,  in  preparation),  MRK  3  shows  a  flat  spectrum  with  a 
photoelectric  cutoff,  consistent  with  a  high  column  density,  but 
not  a  Compton  thick  absorber.  On  the  other  hand,  NGC  1365, 
NGC  3281,  and  NGC  7582  are  most  probably  Compton  thick 
sources  and  the  BAT  data  for  NGC  2992  and  MCG-03-34- 
064  has  a  low  S/N  which  makes  it  difficult  to  constrain  the 
spectrum  well.  These  results  suggest  that  the  low  BAT  to  [O  iv] 
ratio  is  a  marker  for  very  high  column  densities  toward  the 
X-ray  source,  in  agreement  with  previous  studies  (Melendez 
et  al.  2008b;  Rigby  et  al.  2009;  La  Massa  et  al.  2009).  These 
results  argue  in  favor  of  the  unified  model  of  active  galaxies, 
where  the  amount  of  X-ray  emission  suppressed,  relative  to  an 
isotropic  indicator  for  the  AGN  power,  is  related  to  the  absorbing 
column  density  toward  the  X-ray  source  (see  also  Bassani  et  al. 
1999;  Heckman  et  al.  2005;  Netzer  et  al.  2006).  Given  the  high 
ionization  potential  for  Ne+4  (~97  eV),  the  mid-infrared  [Ne  v] 
emission  lines  are  claimed  to  be  an  unambiguous  signatures 
of  an  AGN  (e.g.,  Satyapal  et  al.  2007,  2009).  Therefore,  the 


very  tight  correlation  between  the  L  bat /[O  iv]  and  Z,BAr/[Ne  v] 
14.32  /xm  ratio  confirm  that  [O  iv]  is  AGN-drivenin  this  sample. 

In  Figure  5,  we  plot  both  the  flux  and  luminosity  of  [Nen], 
[Nem]  and  [Oiv]  against  those  of  [Nev]  14.32 /xm.  The  very 
tight  correlation  between  [Nev]  and  [Oiv]  suggests  that  both 
of  these  lines  are  produced  by  the  same  physical  process,  i.e., 
photoionization  by  the  AGN  continuum  (see  also,  Melendez 
et  al.  2008a,  2008b).  This  is  also  true  for  the  tight  correlation 
between  [Ne  v]  and  [Ne  hi]  found  in  our  sample.  While  there 
is  some  scatter  for  [Nem],  which  may  also  be  due  to  star 
formation  (Ho  &  Keto  2007),  the  tightness  of  the  correlation 
suggests  that  [Ne  m]  is  primarily  produced  by  the  AGN  in  these 
objects  (see  also  Gorjian  et  al.  2007;  Goulding  &  Alexander 
2009),  in  agreement  with  the  good  correlation  found  between 
the  [Ne  hi]  and  BAT  luminosity.  Furthermore,  the  tightness  of 
these  correlations,  especially  in  flux-flux,  suggests  that  the 
constant  mid-infrared  ratios  observed  are  not  dominated  by 
aperture  effects.  The  extraction  aperture  for  the  [Oiv]  is  bigger 
than  that  for  the  [Ne  m]  and  [Ne  v]  14  //in,  implying  that  the 
emitting  regions  where  these  lines  originate  are  well  within  the 
central  kpc  of  these  sources,  in  agreement  with  our  previous 
photoionization  studies  (e.g.,  Melendez  et  al.  2008a,  2008b). 
We  found  that  none  of  the  correlations  presented,  e.g.,  between 
the  mid-infrared  lines  and  the  BAT  luminosities,  are  driven  by 
distances  effects,  as  shown  by  the  partial  correlation  analysis 
(see  Table  2)  and  in  agreement  with  the  tightness  of  the  flux-flux 
correlations.  Finally,  the  largest  scatter  occurs  for  [Ne  n],  which 
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Figure  5.  Correlation  between  [Neii],[Nem],  [Nev]  14 /mi,  and  [Oiv]  fluxes 
and  luminosities.  Symbols  are  the  same  as  in  Figure  3.  The  dashed  line  represents 
the  linear  regression  for  the  full  sample  (see  Table  4).  The  statistical  analysis 
and  linear  regression  fits  for  these  correlation  are  presented  in  Tables  2  and  4, 
respectively. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


is  consistent  with  the  interpretation  that  this  line  may  contain  a 
significant  contribution  from  stellar  processes,  resulting  in  the 
enhancement  of  this  emission  line,  as  also  evident  in  Figure  3. 
As  we  discussed  before,  the  [Nev]  emission  lines  are  good 
indicators  for  the  AGN  power,  therefore,  the  good  correlation 
found  between  [Neil]  and  [Nev]  (see  Table  2)  suggest  that,  on 
average,  the  observed  [Ne  n]  emission  in  the  BAT  sample  is  still 
dominated  by  the  AGNs. 

As  we  mentioned  before,  the  ratios  of  high-ionization  lines 
to  low-ionization  lines  can  reveal  the  relative  contributions 
of  the  AGNs  and  star  formation.  In  Figure  6,  we  show  the 
ratios  of  [Nem]/[Neii]  and  [Ne  v]  14/zm/[Nen]  compared 
to  [Oiv]/[Nen]  for  the  BAT  sample.  One  result,  which  was 
also  discussed  in  Melendez  et  al.  (2008a),  is  that  ~62%  and 
~72%  of  our  AGNs  have  [Nem]/[Nen]  and  [Oiv]/[Nen] 
ratios,  respectively,  greater  than  unity.  This  suggests,  again, 
that  in  the  BAT  targets,  the  mid-infrared  emission  lines  are 
dominated  by  the  AGNs.  These  results  are  also  in  agreement 
with  the  high-resolution  Spitzer  analysis  of  the  12  p  m  sample  by 
Tommasin  et  al.  (2008,  2010)  where  they  found  the  same  range 
for  the  [Nem]/[Nen]  and  [Oiv]/[Nen]  ratios.  However,  the 
K-S  test  for  the  [Ne  m]/[Ne  ii]  ratio  between  the  12  /im  sample 
and  the  BAT  sample  returns  a  ~1.9%  probability  of  the  null 
hypothesis,  meaning  that  these  samples  are  statistically  different 
in  their  stellar  content,  with  galaxies  in  the  12  pm  sample  having 
slightly  smaller  ratios  on  average,  1.1  ±  0.9,  than  that  found  for 
the  BAT  sample,  1.3  ±  0.8.  On  the  other  hand,  the  K-S  test  for 
the  [O  iv]/[Ne  in]  ratio  between  the  12  pm  sample  and  the  BAT 
sample  returns  a  ~26.5%  probability  of  the  null  hypothesis. 
These  results  suggest  that  even  the  12  pm  sample  could  be 


log  [O  IV]/[Ne  II] 

Figure  6.  [Ne  m]/[Ne  ii]  and  [Ne  v]14  gm/[Ne  ii]  vs.  [Oiv]/[Nen]  ratios  for 
Seyfert  1,  Seyfert  2,  BAT-detected  AGNs,  and  LINERs  in  the  BAT  sample. 
Symbols  are  the  same  as  in  Figure  3.  The  ~62%  and  ~72%  of  the  objects  in 
our  sample  have  [Ne  m]/[Ne  ii]  and  [O  iv]/[Ne  ii]  have  ratios  bigger  than  unity, 
which  has  been  shown  to  be  typical  of  AGN-ionized  emission-line  gas  (e.g., 
Melendez  et  al.  2008a).  The  [Nev]  and  [Oiv]  ratios  appear  well  correlated,  as 
expected  given  the  tight  correlation  seen  in  Figure  5. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


subject  to  recent  star  formation  contamination,  as  mapped 
through  the  [Nen]  enhancement,  in  agreement  with  previous 
considerations  (Hunt  &  Malkan  1999).  This  mild  contamination 
due  to  star  formation  becomes  noticeable  when  comparing 
the  12  pm  sample  with  the  14—195  keV  sample,  the  latter  of 
which  is  less  biased  toward  star-forming  systems.  Nonetheless, 
the  similarity  on  their  [Oiv]/[Nem]  ratios  suggests  that  both 
samples  have  a  strong  AGN  contribution  to  their  observed 
narrow  line  emission. 

The  tight  correlation  between  [Oiv]  and  [Nev],  evident  in 
Figure  5,  is  seen  here  as  well.  The  [Ne  v]/[Ne  ii]  versus  [O  iv]/ 
[Ne  ii]  plot  shows  a  number  of  sources,  including  several  of  the 
newly  BAT-detected  AGN  and  the  two  LINERs  in  the  sample, 
with  relatively  weak  high-ionization  lines,  which  can  result 
from  a  weak  AGN  (e.g.,  Netzer  et  al.  2004),  strong  nuclear 
star  formation  (e.g.,  Melendez  et  al.  2008a),  extinction  toward 
the  NLR  region  (e.g.,  Netzer  et  al.  2006),  or  shielding  of  the 
NLR  from  the  ionizing  continuum  (e.g.,  Armus  et  al.  2007). 
Comparing  the  [Nem]/[Nen]  to  the  [Nev]  14.32  /xm/[Nen] 
we  found  the  same  overall  trend,  however,  the  slope  of  the  former 
flattens  toward  lower  ratios  of  [O  iv]/[Ne  ii]  possible  due  to  the 
contribution  of  star  formation  to  the  [Nem],  Interestingly,  all 
these  ratios  are  log  normal  distributed  with  a  standard  deviation 
of  only  ~0.3  dex,  except  for  the  [NeV]  14.32  pm/[Q  iv]  with  a 
standard  deviation  of  ~0.2  dex.  The  [O  iv]/[Ne  ii]  emission  line 
ratio  was  also  compared  with  the  redshift  (z)  in  order  to  check 
if  our  previous  results  are  being  affected  by  aperture  affects, 
in  other  words,  if  our  results  are  biased  toward  small  [Nem]/ 
[Nen]  and  [Oiv]/[Nen]  ratios  at  higher  redshifts.  We  found 
that  the  Spearman  rank,  ps  =  0.129;  Pp  =  0.30,  and  Kendall 
test,  r  =  0.08;  /\  =  0.27,  did  not  show  any  correlation  with  z. 
A  similar  result  was  obtained  for  the  [Oiv]/[Nem]  ratios  and 
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Figure  7.  Correlation  between  the  [Ne  v]  24.32  /xm  and  [Ne  v]  14.32  /im  fluxes  and  luminosities  for  our  X-ray  selected  sample.  Symbols  are  the  same  as  in  Figure  3. 
The  low-density  limit  for  this  ratio  is  given  by  the  dashed  line;  AGNs  lying  below  and  to  the  right  have  [Nev]  14.32  /xm/[Ne  v]  24.32  /xm  ratios  smaller  than  the 
theoretical  limit.  The  majority  of  the  newly  detected  BAT  AGNs  have  ratios  below  the  limit,  implying  dust  extinction  toward  the  [Ne  v]14  /xm  emitting  region. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  2 

Statistical  Analysis  for  the  Different  Relationships  Between  the  Mid-infrared  Emission  Lines  and  BAT  Luminosity  for  the  Sample 


Variables 

Ps 

PP 

r 

P 

r 

TP 

<7 

PAP) 

All  Sample 

BAT-[Ne  ii] 

0.66 

<  1  X 

io-6 

0.49 

<  1  X 

IO”6 

0.37 

0.07 

<  1  x  10“6 

BAT-[Ne  hi] 

0.76 

<  1  X 

io-6 

0.58 

<  1  X 

nr6 

0.45 

0.06 

<  1  x  10“6 

BAT-[0  iv] 

0.65 

<  1  X 

ur6 

0.48 

<  1  X 

IO”6 

0.36 

0.06 

<  1  x  IO"6 

BAT-[Ne  v]i4Mm 

0.63 

<  1  X 

ur6 

0.50 

<  1  X 

IO-6 

0.38 

0.07 

<  1  x  IO”6 

[Nem]-[Nen] 

0.89 

<  1  X 

io-6 

0.74 

<  1  X 

1()-6 

0.68 

0.06 

<  1  x  IO-6 

[Oiv]~[Neii] 

0.80 

<  1  X 

io-6 

0.60 

<  1  X 

IO”6 

0.54 

0.06 

<  1  x  IO"6 

[OivHNem] 

0.94 

<  1  X 

ur6 

0.80 

<  1  X 

1()-6 

0.76 

0.06 

<  1  x  IO”6 

[Ne  v]  i4  ^m-PSfe  n] 

0.77 

<  1  X 

ur6 

0.59 

<  1  X 

IO"6 

0.53 

0.06 

<  1  x  10“6 

[Nev]i4/xm-[Neiii] 

0.91 

<  1  X 

io-6 

0.79 

<  1  X 

nr6 

0.75 

0.06 

<  1  x  IO-6 

[Ne  v]i4(im-[Oiv] 

0.94 

<  1  X 

ur6 

0.80 

<  1  X 

IO”6 

0.77 

0.06 

<  1  x  IO-6 

[Nev]24Alm-[Nev]i4 

0.93 

<  1  X 

io-6 

0.80 

<  1  X 

IO"6 

0.77 

0.06 

<  1  x  10~6 

Seyfert  1  Galaxies 

BAT-[Ne  ii] 

0.75 

<  1  X 

io-6 

0.60 

<  1  X 

IO”6 

0.37 

0.11 

<  9  x  IO"4 

BAT-[Ne  hi] 

0.82 

<  1  X 

io-6 

0.66 

<  1  X 

IO”6 

0.41 

0.12 

<  1  x  10-3 

BAT-[Oiv] 

0.72 

<  1  X 

ur6 

0.56 

<  1  X 

io-6 

0.33 

0.11 

<  3  x  10-3 

BAT-[Ne  v] 

0.61 

<  2  x 

ur4 

0.52 

<  1  X 

1()-6 

0.32 

0.12 

<  8  x  10-3 

[OivHNen] 

0.84 

<  1  X 

io-6 

0.66 

<  1  X 

1()-6 

0.54 

0.09 

<  1  x  10“6 

[Oiv]~[Nem] 

0.93 

<  1  X 

io-6 

0.79 

<  1  X 

1()-6 

0.71 

0.08 

<  1  x  10“6 

[Ne  v]  14  ^m-PSfe  n] 

0.77 

<  1  X 

ur6 

0.61 

<  1  X 

IO”6 

0.49 

0.08 

<  1  x  10“6 

[Ne  v]  14  ^nMNe  hi] 

0.86 

<  1  X 

ur6 

0.75 

<  1  X 

1()-6 

0.67 

0.08 

<  1  x  IO”6 

[Ne  v]i4(im-[Oiv] 

0.92 

<  1  X 

io-6 

0.78 

<  1  X 

IO-6 

0.72 

0.08 

<  1  x  10“6 

[Nev]24^m-[Nev]i4 

0.92 

<  1  X 

ur6 

0.81 

<  1  X 

1()-6 

0.77 

0.08 

<  1  x  IO-6 

Seyfert  2  Galaxies 

BAT-[Ne  ii] 

0.63 

<  7.13  x 

ur5 

0.48 

A 

-J 

X 

10  5 

0.40 

0.11 

<  5  x  IO-4 

BAT-[Ne  in] 

0.71 

<  1  X 

ur6 

0.55 

<  1  X 

1()-6 

0.45 

0.10 

<  2  x  10-5 

BAT-[0  iv] 

0.56 

<  6  x 

10~4 

0.42 

A 

X 

1()-4 

0.33 

0.10 

<  2  x  10-3 

BAT-[Ne  v] 

0.65 

<  2  x 

IO”4 

0.49 

<  1  X 

1()-6 

0.38 

0.11 

<  8  x  IO"4 

[Oiv]-[Nen] 

0.76 

<  1  X 

10-6 

0.59 

<  1  X 

1()-6 

0.54 

0.09 

<  1  x  IO-6 

[Oiv]-[Nem] 

0.94 

<  1  X 

ur6 

0.81 

<  1  X 

IO”6 

0.79 

0.08 

<  1  x  10“6 

[Nev]i4Aim-[Nen] 

0.73 

A 

X 

ur5 

0.58 

<  1  X 

1()-6 

0.52 

0.09 

<  1  x  IO”6 

[Nev]i4Aim-[Nein] 

0.94 

<  1  X 

ur6 

0.83 

<  1  X 

1()-6 

0.80 

0.10 

<  1  x  10“6 

[Ne  v]i4(im-[Oiv] 

0.94 

<  1  X 

ur6 

0.80 

<  1  X 

1()-6 

0.78 

0.09 

<  1  x  IO-6 

[Nev]24^m-[Nev]i4 

0.92 

<  1  X 

io-6 

0.80 

<  1  X 

1()-6 

0.78 

0.09 

<  1  x  IO-6 

Notes.  ps  is  the  Spearman  rank  order  correlation  coefficient  with  its  associated  null  probability,  Pp .  r  represents  the  generalized  Kendall’s 
correlation  coefficient  for  censored  data  and  zp  is  the  Kendall’s  coefficient  for  partial  correlation  with  censored  data.  PT  and  Pr(p)  are 
the  null  probabilities  for  the  generalized  and  partial  Kendall’s  correlation  test,  respectively.  We  also  present  the  calculated  variance,  <r,  for 
Kendall  zp.  We  have  used  a  partial  correlation  test  to  exclude  the  effect  of  redshift  (distance)  in  the  luminosity-luminosity  correlations. 


3.2.  Extinction  in  the  Mid-infrared 


the  redshift  where  we  found  no  correlation,  with  a  Spearman 
rank,  ps  =  —0.05;  Pp  =  0.66,  and  Kendall  test,  r  =  —0.05; 
PT  =  0.51. 


If  wavelength-dependent  extinction  is  important  in  the  mid- 
infrared,  it  should  be  evident  when  comparing  the  strengths 
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Figure  8.  [Ne  v]24^m  and  [Ne  ratios  as  a  function  of  host  galaxy  inclination  for  the  BAT  sample.  Symbols  are  the  same  as  in  Figure  3.  The  horizontal  line 

shows  the  low  density  limit.  Note  that  five  out  of  six  of  the  newly  detected  BAT  AGNs  reside  in  inclined  hosts  b/a  <  0.4;  hence,  host  galaxy  extinction  may  be 
responsible  for  their  lack  of  AGN  signatures  in  the  optical. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  3 

Statistical  Analysis  Between  Seyfert  1  and  Seyfert  2  Galaxies 


Parameter 

Seyfert  1 

Seyfert  2 

Pk-s 

(%) 

Measurements 

Available 

Median 

Standard 

Deviation 

Measurements 

Available 

Median 

Standard 

Deviation 

LogBAT/[0  iv] 

38 

2.5 

0.5 

33 

2.4 

0.5 

2.4 

LogBAT/[Nem] 

38 

2.8 

0.4 

33 

2.5 

0.4 

9.0 

LogBAT  /  [Ne  v]  14  32 

32 

3.0 

0.4 

33 

2.9 

0.5 

26.1 

LogBAT/[Nev]24.32,rm 

31 

3.0 

0.5 

30 

2.8 

0.3 

9.5 

LogBAT/ [Nen] 

38 

2.9 

0.4 

33 

2.7 

0.5 

2.8 

[Ne  m]/[Ne  n] 

38 

1.2 

0.1 

33 

1.4 

0.1 

54.8 

[Oiv]/[Nen] 

38 

1.8 

0.4 

33 

1.7 

0.4 

84.5 

[Oiv]/[Nem] 

38 

1.5 

0.1 

33 

1.7 

0.1 

91.4 

[Ne  V]  14.32  /im/[Nev]24.32  fim 

29 

1.0 

0.3 

30 

1.0 

0.4 

94.0 

L  BAT 

38 

43.55 

0.573 

33 

43.28 

0.657 

6.3 

fqNe  ii] 

38 

40.79 

0.613 

33 

40.76 

0.799 

82.6 

^[Ne  v]  14.32  /xm 

32 

40.54 

0.659 

33 

40.33 

0.886 

69.8 

^  [Ne  in] 

38 

40.92 

0.668 

33 

40.74 

0.776 

50.9 

£[Ne  v]24.32  /xm 

31 

40.52 

0.623 

30 

40.44 

0.746 

28.6 

Go  iv] 

38 

40.98 

0.744 

33 

41.00 

0.793 

98.1 

Notes.  The  last  column,  Pk-s>  represents  the  K-S  test  null  probability.  Upper  limits  for  the  [Ne  v]  fluxes  are  not  included.  This  table  also 
includes  information  about  the  numbers  of  Seyfert  1  and  Seyfert  2  galaxies,  median  values,  and  standard  deviations  of  the  mean  for  the 
measured  quantities. 


of  [Nev]  14.32  /xm  and  [Nev]  24.32  /xm.  Dudik  et  al.  (2007) 
calculated  the  theoretical  values  of  this  ratio  as  a  function 
of  density  and  electron  temperature  (see  their  Figure  1).  As 
density  increases,  the  24.32  /xm  line  is  suppressed  relative  to 
the  14.32  /xm  line,  hence,  the  most  likely  way  in  which  the  ratio 
can  fall  below  the  theoretical  low  density  limit  is  extinction, 
suggesting  that  it  should  be  stronger  for  the  shorter-wavelength 
line.  In  Figure  7,  we  show  the  correlation  between  the  two 
[Ne  v]  lines  for  the  BAT  sample,  with  the  theoretical  low  density 
limit  over-plotted.  We  note  that  the  correlation  between  the 
[Ne  v]  emission  lines  is  strong  and  is  similar  in  both  flux  and 
luminosities,  meaning  that  there  is  a  true  linear  correlation  (see 
Table  4). 


For  the  BAT  sample  the  average  [Ne  v]14.32/,m/[Ne  v]24.32|am 
ratio  is  above  the  low  density  limit  and  similar  for  Seyfert  1 
and  Seyfert  2  galaxies,  1.0  ±  0.3  and  1.0  ±  0.4,  respectively. 
The  K-S  test  for  this  ratio  returns  a  ~94.0%  probability  of  the 
null  hypothesis,  i.e.,  that  there  is  no  difference  between  Seyfert 
1  and  Seyfert  2  galaxies;  however,  Seyfert  2  galaxies  show  a 
bigger  dispersion.  For  a  gas  temperature  of  104  K,  these  results 
suggest  an  electron  density  of  ne  ~  103  cm  ’,  in  agreement 
with  previous  studies  (e.g.,  Sturm  et  al.  2002;  Dudik  et  al.  2007; 
Tommasin  et  al.  2008).  On  the  other  hand,  five  out  of  six  of  the 
newly  BAT-detected  AGNs  have  ratios  below  the  low  density 
limit,  with  an  average  of  0.5  ±  0.2,  indicating  possible  dust 
extinction.  In  Figure  8,  we  compared  the  [Ne  v]  line  ratio  to  the 
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Table  4 

Linear  Regressions  for  the  Mid-infrared  and  14-195  keV  Luminosities 


log  X-\og  Y 

Y  =  aX  +  b 

a 

b 

All  Sample 

BAT-[Ne  ii] 

0.72  ± 

0.11 

9.60 

±  4.85 

BAT-[Ne  m] 

0.88  ± 

0.06 

2.58 

±  2.70 

BAT-[0  iv] 

0.93  ± 

0.08 

0.42 

±  3.32 

BAT-[Ne  v]i4/im 

1.00  ± 

0.10 

-3.28 

±  4.51 

[O  iv]-[Ne  ii] 

0.85  ± 

0.07 

5.83 

±  2.67 

[Oiv]-[Nem] 

0.90  ± 

0.04 

3.96 

±  1.65 

[Ne  m]-[Ne  ii] 

0.94  ± 

0.04 

2.22 

±  1.80 

[Nev]14fim-[Nen] 

0.67  ± 

0.06 

13.46 

±  2.35 

[Nev]14/im-[Nem] 

0.82  ± 

0.04 

7.80 

±  1.43 

[Ne  v]i4fim-[0  iv] 

0.93  ± 

0.03 

3.59 

±  1.20 

[Ne  V]24  //.m- [Ne  V]  14 

0.89  ± 

0.05 

4.52 

±  2.20 

Seyfert  1  Galaxies 

BAT-[Ne  ii] 

0.77  ± 

0.08 

7.08 

±  3.69 

BAT-[Ne  hi] 

0.93  ± 

0.08 

0.43 

±  3.68 

BAT-[0  iv] 

0.93  ± 

0.11 

0.37 

±  5.00 

BAT-[Nev]i4(xm 

0.91  ± 

0.15 

0.79 

±  6.65 

[O  iv]-[Ne  ii] 

0.84  ± 

0.10 

6.21 

±  3.91 

[Oiv]-[Nem] 

0.91  ± 

0.07 

3.57 

±  2.78 

[Ne  m]-[Ne  n] 

0.92  ± 

0.06 

3.03 

±  2.56 

[Nev]i4/xm-[Nen] 

0.74  ± 

0.09 

10.65 

±  3.76 

[Nev]14fim-[Nem] 

0.90  ± 

0.06 

4.44 

±  2.53 

[Ne  v]14fim-[0  iv] 

1.00  ± 

0.05 

0.62 

±  2.16 

[Ne  V]24  /xm  [Ne  V]  14 

1.02  ± 

0.05 

-0.69 

±  1.50 

Seyfert  2  Galaxies 

BAT-[Ne  ii] 

0.94  ± 

0.15 

0.15 

±  6.56 

BAT-[Ne  hi] 

0.99  ± 

0.10 

-1.82 

±  4.42 

BAT-[0  iv] 

1.00  ± 

0.13 

-2.17 

±  5.72 

BAT-[Nev]i4(xm 

1.10  ± 

0.12 

-7.10 

±  5.28 

[O  iv]-[Ne  ii] 

0.98  ± 

0.10 

0.48 

±  3.89 

[O  iv]-[Ne  in] 

0.96  ± 

0.05 

1.47 

±  2.14 

[Ne  m]-[Ne  ii] 

1.02  ± 

0.05 

-1.01 

±  2.18 

[Ne  v]14/im-[Nen] 

0.89  ± 

0.09 

4.78 

±  3.62 

[Nev]14/im-[Nem] 

0.87  ± 

0.05 

5.73 

±  2.11 

[Ne  v]14fim-[0  iv] 

0.90  ± 

0.04 

4.45 

±  1.50 

[Ne  v]24^lm— [Ne  v]i4^,m 

1.07  ± 

0.04 

-2.97 

±  1.48 

Notes,  log  X  and  log  Y  represent  the  independent  and  dependent  variables, 
respectively,  a  and  b  represent  the  regression  coefficient  (slope)  and  regression 
constant  (intercept),  respectively.  For  the  relationship  between  the  14-195  keV 
and  mid-infrared  luminosities  we  used  the  ordinary  least-square  regression 
of  the  dependent  variable,  Y,  against  the  independent  variable  X,  OLS(F|X). 
For  the  relationship  between  the  mid-infrared  emission  line  luminosities  we 
used  the  OLS  bisector  method  which  treat  the  variables  symmetrically  (see 
Isobe  et  al.  1990,  for  a  review).  When  censored  data  were  present  we  used 
the  EM  and  the  Buckley-James  (BJ)  method,  if  censored  data  were  present  in 
both  variables  (e.g.,  the  [Nev]i4  —  [Ne  v]24  relationship)  the  Schmitt’s  binned 
method  was  used  (Isobe  et  al.  1986;  Isobe  &  Feigelson  1990).  The  results  from 
the  EM  and  BJ  methods  agree  within  their  estimated  errors;  therefore,  for  the 
sake  of  simplicity,  we  present  only  the  values  from  the  EM  regression. 

host  galaxy  inclination.15  All  of  the  newly  BAT-detected  AGNs 
are  found  in  inclined  host  galaxies  which  argue  in  favor  of  a 
scenario  where  the  optical  signatures  of  the  AGNs  may  have 
been  missed  in  this  group,  in  part,  due  to  extinction  in  the  plane 
of  the  host  galaxy;  as  suggested  by  Kirhakos  &  Steiner  (1990) 
and  Simcoe  et  al.  (1997)  and  more  recently  by  GA09.  However, 
one  must  note  from  Figure  8  that  there  is  no  correlation  between 
the  [Ne  v]  line  ratio  and  the  value  of  the  host  galaxy  inclination. 


15  The  values  for  the  major  and  minor  diameters  of  the  host  galaxy,  a  and  b 
respectively,  were  taken  from  NED. 


Hence,  host  galaxy  inclination  cannot  be  solely  responsible  for 
the  observed  mid-infrared  extinction  in  our  sample. 

In  addition  to  providing  some  insight  into  the  nature  of  the 
newly  BAT-detected  AGNs,  these  plots  show  that  ~34%  of 
Seyfert  1  galaxies  and  ~37%  of  Seyfert  2  galaxies  have  [Ne  v] 
ratios  below  the  low  density  limit,  in  good  agreement  with  the 
percentage  of  Seyfert  galaxies  that  fall  below  the  low  density 
limit  for  the  neon  ratio  from  the  high-resolution  spectroscopy  of 
the  12  jtm  sample  (Tommasin  et  al.  2010).  If  extinction  at  shorter 
wavelengths  were  responsible  for  the  neon  ratio  below  the  low 
density  limit,  it  is  plausible  that  the  tightness  of  the  [Ne  iii]-Lbat 
correlation  may  be  due  to  the  combined  effects  of  mid-infrared 
extinction  and  X-ray  decrement  due  to  Compton  scatter  in 
the  BAT  band  (e.g.,  Melendez  et  al.  2008b).  Similarly,  the 
larger  scatter  for  [O  iv]-Lbat  may  be  due,  in  part,  to  the  lower 
extinction  for  [O  iv],  compared  to  [Ne  hi]  ;  hence,  the  variations 
in  the  observed  X-ray,  due  to  different  column  densities  toward 
the  nuclear  region,  become  the  dominant  effect.  One  must 
note  that  the  mid-infrared  extinction  curved  by  Weingartner  & 
Draine  (2001)  and  Li  &  Draine  (2001;  e.g,  see  their  Figure  16), 
fails  to  predict  the  required  dust  extinction  at  Nn  and 

A24.32iJ.m/NH  to  explain  the  observed  [Ne  v]  ratios  below  the  low 
density  limit.  Other  possibilities  to  explain  the  observed  neon 
ratios  below  the  low  density  limit  must  be  consider,  such  as, 
aperture  effects  and/or  adopted  atomic  data  for  neon;  however, 
these  possibilities  have  also  failed  to  fully  explain  the  observed 
ratios  below  the  low  density  limit  (e.g.,  Dudik  et  al.  2007; 
Tommasin  et  al.  2010). 

4.  COMPARISON  BETWEEN  BAT  AGNs  AND 
OTHER  SOURCES 

Following  the  mid-infrared  properties  for  the  BAT  AGN 
sample  we  investigated  the  relationship  between  different  types 
of  galaxies.  For  comparison  we  searched  the  literature  for  the 
most  recent  Spitzer  high-resolution  IRS  spectroscopy,  from  this 
we  used  the  volume-limited  sample  of  bolometrically  luminous 
galaxies  (LjB  7 >  3  x  109io)  to  a  distance  of  D  <  15  Mpc 
of  GA09  which  includes  Hii  galaxies,  LINERS,  optically 
classified  AGNs  and  optically  unclassified  galaxies.  We  also 
used  the  recent  atlas  of  starburst  galaxies  by  B09  that  includes 
starburst  galaxies  with  and  without  evidence  of  an  AGN.  Finally, 
we  used  a  sample  of  extreme  starburst  galaxies  which  includes 
blue  compact  dwarfs  (BCD)  from  H09.  When  a  source  was 
present  in  more  than  one  compilation  we  gave  priority  to  the 
data  analysis  and  extraction  presented  in  GA09,  followed  by 
B09  and  H09. 

Using  [Ne  v]  14.32  fim,  GA09  identified  AGN  activity  in 
sources  that  were  previously  not  classified  by  their  optical 
spectra  as  AGNs.  While  theoretical  models  predict  that  some 
fraction  of  [Ne  v]  can  be  produce  by  energetic  starbursts  (e.g., 
Abel  &  Satyapal  2008),  H09  did  not  detect  [Nev]  emission 
in  their  sample  of  extreme  starburst  galaxies.  However,  some 
problems  arise  when  looking  at  the  [Nev]  emission  in  AGNs. 
As  we  mentioned  before,  [Nev]  14.32 /rm  and  24.32 pm 
were  not  detected  in  all  the  sources  in  our  sample  of  BAT 
AGNs.  This  suggests  that  because  of  dust  extinction  and 
instrumental  sensitivity,  [Nev]  emission  could  go  undetected 
even  in  intrinsically  luminous  AGNs.  On  the  other  hand,  [Ne  ii], 
[Nem],  and  [Oiv]  are  present  in  AGNs,  BCD,  starburst,  Hu 
galaxies  and  have  a  wide  range  of  ionization  potentials  and 
critical  densities  which  allow  us  to  study  the  connection  between 
the  active  nucleus  and  star  formation.  One  must  note  that  H09 
found  that  the  most  likely  contributor  to  the  [Oiv]  emission 
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Figure  9.  [Nem]/[Neii]  vs.  [Oiv]/[Nem]  for  a  sample  of  different  sources,  including  our  BAT  AGNs,  LINERs,  starburst  galaxies,  Hu  galaxies,  and  blue  compact 
dwarf  from  high-resolution  Spitzer  spectra  from  Goulding  &  Alexander  (2009,  GA09),  Bemard-Salas  et  al.  (2009,  B09),  and  Hao  et  al.  (2009,  H09).  Note  that  symbols 
are  different  than  in  previous  figures,  see  legend.  The  dashed  lines  represent  a  value  of  unity  for  the  [Ne  m]/[Ne  n]  and  [O  iv]/[Ne  hi]  ratios.  The  rectangles  represent 
the  boundaries  for  the  sample  of  74  ULIRGs  (solid  line)  and  34  PG  quasars  (point-dashed  line)  from  Veilleux  et  al.  (2009).  Interestingly,  half  of  the  BAT  sample  can 
be  uniquely  distinguished  from  SB /H  n/BCD  galaxies  by  having  both  [Ne  m]/[Ne  n]  and  [O  iv]/[Ne  hi]  ratios  greater  than  unity.  For  the  ratio  errors,  we  propagated 
the  uncertainty  for  individual  emission  lines  in  quadrature.  The  error  bars  at  the  bottom  right  comer  show  an  average  error  for  the  ratios  presented  in  this  plot. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


in  starburst  and  BCD  galaxies  is  Wolf-Rayet  stars.  In  this 
regard,  and  despite  the  fact  that  [Oiv]  emission  traces  the 
AGN  intrinsic  luminosity  in  pre-selected  AGNs  (Melendez  et  al. 
2008b),  [O  iv]  cannot  be  solely  associated  with  the  power  of  the 
AGNs. 

In  Figure  9,  we  compare  the  [Neni]/[Nen]  versus  [Oiv]/ 
[Ne  hi]  ratios.  From  this  comparison  there  are  four  distinctive 
regions  to  note:  (1)  the  BAT  AGNs  branch;  (2)  below  the  BAT 
AGNs  are  the  starburst  (SB)  and  Hii  galaxies  which  have  an 
apparent  [Neil]  excess;  (3)  above  and  to  the  left  of  the  BAT 
AGNs  are  the  BCD  and  two  WR  galaxies  from  GA09  (II  Zw40 
and  NGC  1569),  and  (4)  LINERs,  which  seem  to  follow  a 
connecting  path  between  AGNs  and  SB/Hn  galaxies  but  are 
mainly  found  in  the  lower-luminosity  region  of  the  BAT  AGNs 
branch.  Despite  the  fact  that  [Oiv]  has  been  detected  in  BCD, 
SB,  and  H  ii  galaxies,  WR  and  BCD  galaxies  cannot  produce 
[Oiv]  as  effectively  as  [Nem],  unlike  in  AGNs,  making  these 
sources  have  relative  low  [Oiv]/[Nem]  ratios.  On  the  other 
hand,  the  hardness  of  the  photoionizing  continuum  from  WR 
stars  and  other  energetics  phenomena,  such  as  Ultraluminous 
X-ray  Sources  (ULXs),  can  ionize  Ne+  into  Ne2+  resulting  in 
relatively  weak  [Nen]  (e.g.,  Berghea  et  al.  2010),  thus  the  high 
[Ne  m]/[Ne  ii]  ratios  observed  in  some  BCD/SB  galaxies. 

There  is  also  a  small  “merge”  region,  at  the  lower  ionization 
end  of  the  BAT  AGNs  branch,  that  comprises  LINERs  and 
starburst  galaxies  with  evidence  for  an  AGN.  This  region 
exemplify  the  different  properties  of  LINERs,  e.g.,  IR-luminous 


and  faint  LINERs,  with  the  latter  having  similar  SED  than  that 
found  in  starburst  galaxies  (Sturm  et  al.  2006).  Moreover,  in 
this  region  we  found  five  of  the  [Ne  v]  AGNs  from  GA09,  some 
of  which  do  not  have  optical  signatures  for  an  AGN.  Among 
these  [Ne  v]  AGNs,  NGC  660  and  NGC  3628  have  shown  some 
evidence  to  harbor  a  LINER  (Filho  et  al.  2004;  Dudik  et  al.  2005; 
Flohic  et  al.  2006).  Also,  based  on  this  diagnostic,  we  found  that 
NGC  520,  NGC  4676,  and  NGC  1614  could  harbor  an  AGN, 
although  these  starburst  galaxies  do  not  show  optical  evidence 
of  such  (B09).  There  is  X-ray  evidence  for  a  Compton  thick 
AGN  in  NGC  1614  (Risaliti  et  al.  2000a),  in  agreement  with 
our  diagnostic.  Interestingly,  these  possible  AGNs  are  found 
in  interacting  systems  or  mergers  with  strong  star  formation. 
Perhaps,  interacting  or  merger  systems  may  provide  enough 
fuel  to  the  nuclear  regions  of  the  galaxy  to  trigger  both  the 
nuclear  star  formation  and  the  AGNs. 

We  also  over  plot  in  Figure  9  the  sample  of  74  ULIRGs 
and  34  Palomar-Green  quasars  (PG  QSOs),  observed  with  the 
IRS  Spitzer ,  from  Veilleux  et  al.  (2009).  From  this  comparison 
it  is  clear  that  ULIRGs  show  a  wide  range  of  AGNs  and  SB 
contribution  to  their  mid-infrared  emission  lines,  filling  the  gap 
between  the  BAT  AGNs  and  the  SB /H  ii  branch  and  overlapping 
the  small  “merge”  region  defined  above,  in  agreement  with  the 
idea  that  ULIRGs  are  composite  systems  mainly  powered  by 
stellar  activity  (Armus  et  al.  2007;  Veilleux  et  al.  2009).  On 
the  other  hand,  PG  QSOs  overlap  with  the  BAT  AGN  branch, 
in  agreement  with  the  high,  typically  larger  than  ~80%,  AGN 
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contribution  to  their  bolometric  luminosities  (Veilleux  et  al. 
2009).  One  must  note  that  there  are  no  ULIRGs  detected  in  the 
Swift-BAT  catalog  (Tueller  et  al.  2010;  Cusumano  et  al.  2010), 
except  for  the  composite  system  NGC  6240  (Armus  et  al.  2006). 

Moreover,  most  of  the  non-BAT  AGNs  presented  in  this  work, 
e.g.,  the  SB+AGN  sources  from  B09  and  [Nev]  14.32 /rm- 
AGN  (GA09),  in  other  words,  AGNs  that  have  been  selected 
because  their  optical  and  mid-infrared  emission  line  properties 
have  smaller  mid-infrared  ratios  than  that  found  for  the  BAT 
AGNs.  In  fact,  half  of  the  BAT  sample  can  be  uniquely  distin¬ 
guished  from  SB/Hn/BCD  galaxies  by  having  both  [Nem]/ 
[Nen]  and  [O  iv]/[Ne  hi]  ratios  greater  than  unity.  This  result 
suggests  that  the  BAT  sample  represents  a  unique  opportu¬ 
nity  to  study  high-ionization  AGNs,  i.e.,  sources  in  which  their 
optical/mid-infrared  emission  signatures  are  dominated  by  the 
AGNs.  Figure  9  shows  that  the  [Ne  m]/[Ne  n]  ratio  as  a  stand¬ 
alone  diagnostic  may  overestimate  the  fraction  of  AGNs  be¬ 
cause  Wolf-Rayet  stars  or  another  energetic  phenomena  could 
produce  [Nem]/[Neii]  ratios  even  higher  than  those  found  on 
the  BAT  AGNs.  However,  by  further  constraining  this  ratio  by 
using  [Oiv]  we  were  able  to  separate  AGNs  from  energetic 
starbursts  and  other  sources. 

As  we  mentioned  before,  caution  must  be  taken  when  com¬ 
paring  fluxes  between  Spitzer  high-resolution  orders  because  of 
the  differences  in  aperture  size.  However,  given  the  similarities 
between  the  ionization  potentials  for  [Nem]  and  [Oiv],  it  is 
likely  that  both  lines  originate  in  similar  regions,  thus  intro¬ 
ducing  less  uncertainty  in  the  use  of  the  [Oiv] /[Nem]  ratio. 
Nevertheless,  in  Hii/SB/BCD  galaxies  we  could  be  underesti¬ 
mating  the  amount  of  [Ne  hi]  that  is  associated  with  the  same 
physical  conditions  that  produce  [O  iv],  e.g.,  WR  stars.  This  will 
cause  this  ratio,  when  derived  from  equal  apertures  extractions, 
to  be  displaced  toward  even  smaller  values  thus,  enhancing  the 
differences  between  these  sources  and  our  sample  of  hard  X-ray 
selected  AGNs. 

5.  CONCLUSIONS 

Using  high-resolution  Spitzer  IRS  spectra,  we  have  examined 
the  mid-infrared  emission-line  properties  of  a  sample  of  hard 
X-ray  selected  AGNs,  detected  by  Swift /BAT.  Our  principle 
conclusions  are  as  follows. 

1.  The  luminosity  distribution  for  the  mid-infrared  emission 
lines  and  BAT  continuum  luminosities  show  no  differ¬ 
ences  between  Seyfert  1  and  Seyfert  2  populations  for 
the  BAT  sample.  The  correlations  between  all  the  mid- 
infrared  emission  lines  and  BAT  in  both  flux-flux  and 
luminosity-luminosity  are  statistically  significant,  even 
when  factoring  the  distance  effect  in  luminosity-luminosity 
correlations.  The  dispersion/tightness  in  these  correlations 
is  due  to  differences  in  the  X-ray  absorbing  column  densi¬ 
ties,  dust  extinction,  and/or  nuclear  star  formation  activity. 
Moreover,  the  tight  correlation  found  in  the  [Nem]-BAT 
relationship  suggests  that,  on  average,  there  is  no  strong  en¬ 
hancement  due  to  star  formation  in  the  [Ne  hi]  emission  in 
the  BAT  sample.  Also,  the  slopes  for  the  [Ne  m],[Ne  v],  and 
[O  iv]  versus  BAT  luminosities  relationships  are  smaller  in 
Seyfert  1  galaxies  than  in  Seyfert  2s  (which  are  around 
unity),  which  suggests  that,  while  the  amount  of  extinction 
toward  the  NLR  is  similar  in  both  types,  the  X-ray  absorb¬ 
ing  columns  are  large  enough  in  Seyfert  2s  to  affect  the 
hard  X-ray  band,  confirming  the  results  of  Melendez  et  al. 
(2008b)  and  Rigby  et  al.  (2009).  This  result  is  in  agreement 


with  the  fact  that  the  BAT/[0  iv]  ratio  statistically  separates 
Seyfert  Is  and  Seyfert  2s. 

2.  Although  all  of  the  correlations  among  the  mid-infrared 
emission  lines  are  strong,  the  worst  correlations  are  for 
[Nev]-[Nen]  and  [Oiv]-[Nen],  because  of  enhance¬ 
ment  of  the  [Nen]  from  nuclear  stellar  activity  (see  also 
Melendez  et  al.  2008a).  While  the  tightness  of  these  mid- 
infrared  correlations  suggests  that  dust  extinction  is  not  the 
driving  physical  process  behind  the  mid-infrared  relation¬ 
ships,  approximately  ~40%  (including  upper  limits)  of  the 
sample  have  values  for  the  ratio  of  the  [Ne  v]  emission  lines 
below  the  low-density  theoretical  limit,  suggesting  dust  ex¬ 
tinction  as  the  physical  process  responsible.  Exploring  this, 
we  found  that  all  of  the  newly  discovered  BAT  AGNs  in  our 
sample,  which  are  under-luminous  in  [Oiv]  and  [Nev]  14/ 
24  p m,  are  found  on  inclined  host  galaxies,  and  all  but  one 
have  [Ne  v]  ratios  below  the  critical  density  limit.  Hence, 
it  is  likely  that  the  newly  found  BAT  AGNs  in  our  sample 
lack  optical  AGN  signatures  because  of  host  galaxy  extinc¬ 
tion  toward  their  NLRs.  However,  the  lack  of  correlation 
between  host  galaxy  inclination  and  the  neon  ratios  suggest 
that  extinction  along  the  plane  of  the  host  galaxy  cannot  be 
responsible  for  the  observed  extinction  in  all  the  BAT  AGN 
samples. 

3.  We  compared  the  BAT  AGNs  with  different  starburst  and 
H  ii  galaxies,  so-called  [Ne  v]  active  galaxies,  and  LINERs 
(Goulding  &  Alexander  2009;  Bernard-Salas  et  al.  2009; 
Hao  et  al.  2009).  We  found  that  the  BAT  AGNs  fall  into  a 
distinctive  region  based  on  the  [Ne m]/[Ne ii]  and  [Oiv]/ 
[Nem]  ratios.  Using  [Nem]  and  [Oiv]  emission,  previ¬ 
ously  connected  with  AGN  power  (e.g.,  Gorjian  et  al.  2007; 
Melendez  et  al.  2008b),  does  not  unambiguously  identify 
AGNs  as  a  stand-alone  diagnostic  because  Wolf-Rayet  stars 
or  another  energetic  phenomena  (perhaps  ULXs)  could  en¬ 
hance  the  observed  emission.  While  it  is  likely  that  de¬ 
tection  of  [Ne  v]  indicates  the  presence  of  an  AGN,  the 
strongest  of  the  [Ne  v]  lines  have  ~  1  /3  less  flux  than  [O  iv] 
and  thus  will  be  more  difficult  to  detect  in  weak  or  faint 
AGNs.  Therefore,  our  composite  method  using  the  [Ne  ii], 
[Nem],  and  [Oiv]  represents  a  strong  and  simple  diagnos¬ 
tic  by  using  only  three  emission  lines  to  identify  an  AGN. 
Based  on  this,  we  found  that  NGC  520,  NGC  4676,  and 
NGC  1614  could  harbor  an  AGN,  although  these  starburst 
galaxies  do  not  show  optical  evidence  of  such  (B09).  There 
is  X-ray  evidence  for  a  Compton  thick  AGN  in  NGC  1614 
(Risaliti  et  al.  2000a),  in  agreement  with  our  diagnostic. 
Such  line  diagnostic  will  be  particularly  useful  to  analyze 
spectra  from  new  IR  missions,  such  as  the  James  Webb 
Space  Telescope  (Gardner  et  al.  2006). 

We  also  found  that  ULIRGs  and  PG  QSOs  occupy  two  dis¬ 
tinctive  regions  in  our  emission  line  diagnostic.  Most  ULIRGs 
fall  into  the  gap  between  the  BAT  AGNs  and  the  SB /H  ii  branch, 
in  agreement  with  the  idea  that  ULIRGs  are  composite  systems 
mainly  powered  by  stellar  activity  (Armus  et  al.  2007;  Veilleux 
et  al.  2009).  On  the  other  hand,  PG  QSOs  overlap  with  the 
BAT  AGN  branch,  in  agreement  with  the  high,  typically  larger 
than  ~80%,  AGN  contribution  to  their  bolometric  luminosities 
(Veilleux  et  al.  2009). 

Finally,  most  of  the  non-BAT  AGNs  presented  in  our  study, 
AGNs  that  have  been  selected  because  their  optical  and  mid- 
infrared  emission  line  properties,  have  smaller  mid-infrared 
ratios  than  that  found  for  the  BAT  AGNs.  In  this  regard,  half  of 
the  BAT  sample  can  be  uniquely  distinguished  from  SB/Hii/ 
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BCD  galaxies  by  having  both  [Ne  m]/[Ne  n]  and  [O  iv]/[Ne  hi] 
ratios  greater  than  unity.  Moreover,  when  comparing  the  12  /xm 
and  our  BAT  selected  AGNs  we  found  that  the  [Nem]/[Nen] 
ratio  distribution  between  the  samples  is  statistically  different 
with  sources  in  the  12  /xm  sample  having  on  average  lower  ratios 
than  that  found  in  the  BAT  AGNs,  or  alternatively  higher  recent 
stellar  activity.  This  mild  contamination  due  to  star  formation 
becomes  noticeable  when  comparing  the  12  /xm  sample  with  the 
14-195  keV  sample,  the  latter  of  which  is  less  biased  toward 
star-forming  systems.  Despite  the  fact  that  both  samples  have  a 
strong  AGN  contribution  to  their  observed  narrow  line  emission, 
this  result  suggests  that  the  BAT  sample  represents  a  unique 
opportunity  to  study  high-ionization  AGNs,  sources  in  which 
their  optical/mid-infrared  emission  signatures  are  dominated 
by  the  AGNs,  thus,  providing  the  most  representative  sample  in 
terms  of  galaxy  population  and  stellar  content. 
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